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A Widely-Separated, Highly-Occluded Companion to the Nearby Low-Mass T 

Tauri Star TWA 30 1 
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ABSTRACT 

We report the discovery of TWA 30B, a wide (~3400 AU), co-moving M dwarf 
companion to the nearby (~42 pc) young star TWA 30. Companionship is confirmed 
from their statistically consistent proper motions and radial velocities, as well as a 
chance alignment probability of only 0.08%. Like TWA 30A, the spectrum of TWA 
30B shows signatures of an actively accreting disk (H I and alkali line emission) and 
forbidden emission lines tracing outflowing material ([O I], [O II], [O III], [S II], and 
[N II]). We have also detected [C I] emission in the optical data, marking the first such 
detection of this line in a pre-main sequence star. Negligible radial velocity shifts in the 
emission lines relative to the stellar frame of rest (AV < 30 km s -1 ) indicate that the 
outflows are viewed in the plane of the sky and that the corresponding circumstellar disk 
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is viewed edge-on. Indeed, TWA 30B appears to be heavily obscured by its disk, given 
that it is 5 magnitudes fainter than TWA 30A at K-band despite having a slightly earlier 
spectral type (M4 versus M5). The near-infrared spectrum of TWA 30B also evinces 
an excess that varies on day timescales, with colors that follow classical T Tauri tracks 
as opposed to variable reddening (as is the case for TWA 30A). Multi-epoch data show 
this excess to be well-modeled by a black body component with temperatures ranging 
from 630 to 880 K and emitting areas that scale inversely with the temperature. The 
variable excess may arise from disk structure such as a rim or a warp at the inner disk 
edge located at a radial distance of ~ 3-5 Rq. As the second and third closest actively 
accreting and outflowing stars to the Sun (after TWA 3), TWA 30AB presents an ideal 
system for detailed study of star and planetary formation processes at the low-mass end 
of the hydrogen-burning spectrum. 

Subject headings: Galaxy: open clusters and associations: individual (TW Hydrae 
Association) - stars: circumstellar matter, stars: pre-main-sequence - stars: evolution 
- stars: individual (2MASS J11321822-3018316, TWA 30A, TWA 30B) - stars: low- 
mass, brown dwarfs 



Introduction 



The TW Hydrae Association (TWA) comprises the nearest population of classical T Tauri stars 
(cTTS), yet its membership r emains incomp l ete. Formally recognized as a young physical associa- 
tion composed of five st ars by Kastn er et al.1 (119971') . its me mbership swelled to 24 systems in little 



more than a decade (See lMamajekll2005l ; lLooper et al.ll2010l and references t here- in). Many of these 
systems have been intensely studied because of their young age (~8 Myr; e.g., Soderblom et al.ll998l: 
Webb et al. 1999 : Barrado v Navascues 2006 ) and close proximity to the Sun f^ 50 pc: llVlamaiek 
2005). The natal cloud from which these young stars sprung has since dispersed (jTachihara et al. 
2009). Mostly low-mass stars and brown dwarfs (<0.08 M Q ) populate TWA, with only one object 
having a spectral type earlier than K5: the AO V star TWA 11A. TWA members have a variety of 
disk properties at a critical era in protoplanetary formation. Therefore, identification of additional 
members provides invaluable case studies of disks and outflows of low-mass objects in our Galactic 
backyard. 

Associations typically have an overdensity of stars against the background stellar population 
that assist in membership identification. However, the density of TWA members is roughly two 
orders of magnitud e less than the local disk density: ~0.002 M Q pc -3 versus ~0.1 M Q pc~ 3 



(jvan Leeuwenll2007l ) . The loose grouping of TWA members and their proximity to the Sun translate 
into a projected extent of ~1000 square degrees on the sky. Hence, the wide expanse of the TWA 



requires other methods for identifying new members. Past technic 



ray brightness in ROSAT All-Sky Survey data (RASS: IWebb et al- 



lies include identification by X - 



1999 



Zuckerman et alJl2001bl ) 
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IR excess in the IRAS Point Source Catalog (Ide la Reza et al.lll9 89; Grcgori o-Hetem et al.lll992l ) 



2001 



with spectroscopic 



and kinematic se lections from the Tycho-2 catalog (jMakarov &: Fabriciusl ; 
confirmation by ISong et al.l 120020 . These surve ys list only very bright objects, i.e., Tycho-2 is 



statistically complete down to only V=11.5 mag ([Hog et al.l l2000). Therefore, such programs have 
been heavily biased against late-type members of the TWA. 



Nevertheless, there have been several successful searches for substellar TWA members. Gizis 



(|2002l ) conducted a NIR photometric search over the expanse of TWA using the Two Micron All-Sky 
Survey Point Source Catalog (2MASS PSC. ICutri et alJbood : Iskrutskie et alJbood ) and identified 
the two M8 substellar members 2MASS J113 9-3159 (TWA 26) and 2M ASS J1207-3932 (TWA 



27). High resolution imaging o f TWA 27A (IChauvin et al.l 12004 



low-mass companion (~8 M r- T : iMohanty et al.l 120071 ). Previously, 



2005 ) revea l ed a late-L type 



Webb et al.1 ()1999l ) performed 



a high-resolution imaging survey of known TWA members and found the M8 substellar member 
TWA 5B. A proper m otion survey of the Southern sky using the SuperCOSMOS Sky Survey (SSS; 



Hamblv et al1l2001al |bl|3) uncovered the M8 substellar member SSSPM J1102-3159 (|Scholz et al 



200 



5|; TWA 28). Likewise, a near-infrared (NIR ) photometric search for nearby, ultracool dw arfs 



using the DEep Near- Infrared Survey (DENIS; lEpchtein et al.lll997l ) by lLooper et al.l (|20071 ) un- 
covered the M9.5 substellar member DENIS J1245-4429 (TWA 29). 



Looper et al.1 (|2010l ) recently identified another new member of TWA from a survey which 
pairs red candidates in 2MASS PSC data with ROSAT X-ray catalog data - the highly variable 
M5 db 1 pre-main sequence (PMS) star TWA 30. In this paper we report the serendipitous discovery 
of a widely separated companion to TWA 30: 2MASS J11321822-3018316, which we term TWA 
30B. We describe our identification of this source and follow-up observations in Section 2, as well 
as additional spectroscopic coverage of TWA 30A. In Section 3, we review our kinematic and 
spectroscopic analysis of TWA 30B. In Section 4, we propose a toy model for the origins of the NIR 
excess emission in TWA 30B, compare the properties of TWA 30A and B, and discuss the possibility 
that this wide binary is a member of the nearby f3 Pic Association. We give our conclusions in 
Section 5. 



2. Observations 
2.1. Discovery 

We serendipitously discovered 2MASS Jl 1321822— 3018316 while making color finder charts 
for TWA 30A (see Figured]). 2MASS PSC data showed this object as nearby (~80") and red 
(J - K s ~L6). Using sky-survey plates from USNO-B1.0 and 2MASS (Table HJ, we found a 
similar proper motion for this candidate companion to TWA 30A (See Section 3.1). Follow-up 
spectroscopy and astrometric measurements confirm its similar signatures of youth and shared 
kinematic properties. From these multiple lines of evidence of companionship, we term this object 
TWA 30B. 
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2.2. Optical Spectroscopy: Magellan/MagE 



TWA 30B was observed twice with the Magellan Echellete (MagE; Marshall et all 120081 1 on 
the 6.5 m Clay Telescope at Las Campanas Observatory, a cross-dispersed optical spectrograph 
which covers the ~3000-10500 A passband in a single exposure. Conditions on the first night were 
clear, while on the second night thin clouds were present (all observations are summarized in Table 
[2]). The tf!7 slit was employed and aligned to parallactic for each observation, resulting in R = A 
/ AA ~ 4100 with no binning. A flux standard was observed during each run, and ThAr arc 
calibration images were obtain ed after each science ex posure. The spectrum was reduced with the 
MASE reduction tool pipeline (jBochanski et al.1 12009). The reduced spectra from both epochs are 
shown in Figure [2j 



2.3. NIR Spectroscopy: IRTF/SpeX 



TWA 30B was observed over one year with 18 epoc hs of low-resolution near-infrared spec- 
troscopy with the SpeX spectrograph (jRavner et al.ll2003l ) on the 3.0 m NASA Infrared Telescope 
Facility (IRTF). The prism-dispersed mode of SpeX was used in conjunction with the 0'.'5 slit, 
providing 0.7-2.5 ^m spectroscopy at a mean resolution of R ~150. We also observed TWA 30B 
three times with the SpeX spectrograph in cross-dispersed mode (SXD) with the / /5 slit, resulting 
in a mean resolution of R ^1200 over ~0.8-2.4 /im. The details of these observations are recorded 
in Table El Of note, all nights but 2010 January 24 UT and 2010 May 23 UT were clear and 
stable. For telluric correction and flux calibration, observations of a nearby AO V star were made 
immediately after each observation along with internal flat-field and argon arc la mp frames for 



calib r ations. All reduct ions were carried out with the Spextool package version 3.4 (jCushing et al 
2004!: 



Vacca et al.ll2003l ) using standard settings. 



In addition, we observed TWA 30A over eight epochs spanning six months using the SpeX 
spectrograph in SXD mode with the / /5 slit (see Table [2]). This NIR data complements our previous 
spectroscopic coverage of TWA 30A (Looper et al. 2010) by extending the baseline coverage by 11 
months. Observations and reductions were carried out in the same manner as described above for 
TWA 30B. 



2.4. Imaging 



TWA 30A was observed four times from 2008 December to 2 010 January at the Cerro Tololo 
4.0 m Blanco Telescope with the Infrared Side Port Imager (ISPI; Ivan der Bliek et alJl2004T). using 
the J-band filter a s part of an ongoing brown dwarf astrometric program (IFaherty et al.1 |2009j; 



Faherty et al.ll2010l ). While TWA 30 A was the primary target, the large field of view (fov) of the 



detector (~10 arcmin fov with a 0'.'3 pixel plate scale) allowed us to also derive positions for TWA 
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30B (equinox J2000; see Table [T]). Images were taken in an identical manner: three images of 10 s 
exposures using three co-adds at each of the three dither positions. At the beginning of the night, 
we obtained dark frames and domeflats. Raw images were median-combined to produce sky frames, 
which were then subtracted from the raw data. Otherwise, reduction procedures were based on the 
prescriptions put together by the ISPI tea njj], utilizin g a combination of routines from the Image 
Reduction and Analysis Facility (IRAP@; iTodvl 1 19861 ) . as w ell as the pub l icly a vailable software 
packages: WCSTOOLS and SWARP. For further details, see lLooper et all (hoioh . 



3. Analysis 
3.1. Kinematics 



We determined a proper motion for TWA 30B of fj, a = — 83 ± 9 masyr -1 and [is = — 30 ± 9 
masyr -1 , using a weighted mean of the catalog positions of this object spanning over thirty years 
combined with our own measurements (Table [T]) . SSS lists an indepe ndent measurement of jx a = 
—80 ± 11 masyr -1 and fj,g = —35 ± 9 masyr -1 ( Hamblv et alj|2001al ). We adopted our calculated 
proper motion for analysis. These measurements are consistent within their respective errors to the 
proper motion mea surement of TWA 30A, which has fj, a = —89.6 ± 1.3 masyr -1 , fig = —25.8 ± 1.4 
masyr -1 (UCAC3: [Finch et alJl201oh . 



Similarly, the motion of T WA 30B is consistent with TWA membership based on a convergent 
point analysis (jMamajekll2005l ). Assuming a group vector of ( U, V, W)=(— 10.2 ± 0.5, —18.3 ± 0.5, 
—4.9 ± 0.5) for TWA (Mamajek, in prep), the proper motion of TWA 30B towards the convergent 
point is fj, v = 88.3 ± 9.0 masyr -1 and the perpendicular proper motion is \i T = 0.0 ± 9.0 masyr -1 , 
i.e., perfectly consistent with motion toward the convergent point. The predicted distance from this 
analysis is 44.0 ± 4.9 pc from the Sun. At its celestial position and estimated distance, TWA 30B 
resides only 5.4° or 9.7 pc from the centroid for the TWA (a = 174.2°, S = -35.6°, d = 51.5 ± 2.5 
pc; Mamajek, in prep). 

We calculated the radial vel ocities (RV) of both epo chs of the MagE spectra by cross -correlation 

with t he M4 SPSS templa te from lBochanski et al.l ([20071 ) using the xcorl package in IDL ([Mohantv Basri 



2003 



West fc Basri hoOfll ). Our analysis covered the wavelength r 



ange: 



6800-8900 A, a region 

that encompasses the TiO molecular bandhead that begins at 7056 A. We find RVs of 10 ± 4 
km s -1 and 13 ± 4 km s -1 for epochs 2009 March 4 and May 30, respectively, averaging to 12 ± 3 
km s -1 . This is consistent with the predicted RV of a TWA member at the position of TWA 30B 
(11.0 ± 1.0 km s -1 ). 



See http://www.ctio.noao.edu/instruments/ir_instruments/ispi/ 



2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Assuming the two components constitute a physical pair we derive properties of the system 
(see Table [3]). The weighted mean distance of the pair is 42 ± 2 pc. At an angular separation of 
8tf.'2, this translates into a projected physical separation of p = 3370 ± 1 60 AU. Statistically, th e 
actual phy sical separation can be estimated as 1.26 x p« 4250 ± 200 AU ([Fischer Marcvlll992l ). 
Using the iBaraffe et al.l (|1998l ) models and assuming an age of 8 Myr, we estimate the physical 
parameters for TWA 30A (M5) as log T e// = 3.49, M = 0.11 M and for TWA 30B (M4) as log 
T e ff = 3.51, M = 0.19 Mq. Hence the binary's total mass is ~0.30 M . Assuming an actual 
semi-major axis of 4250 AU and a circular orbit, this yields ia = 1560 AU, va = 0.09 km s _1 and 
tb = 2690 AU, vb = 0.16 km s . The period for this configuration is ~0.5 Myr, translating into 
~16 orbits over the assumed age of the system. 

An interesting question arises - is this system really bound or simply a chance alignment of 
two TWA members? Assessing the boundedness for these systems presents difficulties as their pre- 
dicted orbital velocities are on the order of ~0.1 km s , whereas our current astrometric accuracy is 
greater. To calculate the chance alignment of two members, we generated a Monte Carlo simulation 
of 10,000 fake TWA clusters given the following distributions of RA, Dec, and distance with normal 
deviates generated for each parameter: (RA = 174.6°, la dispersion = 11.9°), (Dec = —36.3°, la 
dispersion = 6.3°), (distance = 54.7 pc, la dispersion = 12.5 pc). The mean distance and disper- 
sion are based on revised cluster parallax distances to the individual members and are calculated 
based on the best available proper motions and a new space velocity and convergent point solution 
(RA = 99.8°, Dec = —27.7°, Stot = 22.0 km s _1 ). We find the median minimum separation be- 
tween two TWA members as 2816", with a 95% confidence level of 520-6433" and a 68% confidence 
level of 1423-4548". In only 8 simulations out of 10,000 was the minimum separation between two 
TWA members <80". We therefore estimate a 0.08% probability of a chance alignment for the 
components of TWA 30 AB. 

Of note, TWA contains at least one other proposed very wide binary - TWA 1 and 28, which 
has an angular separation of 735.3" or a projected physical separation of 40,440 AU at a distance 
of 55 pc. The co-motion between TWA 1 and 28 is 1.0 ± 1.9 masyr -1 , which translates into an 
orbital motion of 0.3 ± 0.7 km s . In only 4.74% of our simulations is the minimum separation 
between fake TWA members less than 735.3". The astro metric data for TWA 1/28 provides further 
support f or the pair to b e co-di stant and co-moving. The Ivan Leeuwenl (120071 ) astrometry for TWA 



1 and the lTeixeira et al.1 ( 20091 ) astrometry for TWA 28 agree within 1.0 ± 1.9 masyr 1 in fj, a and 



0.1 ± 1.6 masyr -1 in //^. This translates into a difference of 0.27 ± 0.66 km s -1 if the pair are 
co-distant. If they share velocity but not distance, this places an upper limit on the difference in 
their distance at 3.1 pc (3<r). These simulations suggest (along with precise astrometry for TWA 1 
and 28) that both the TWA 30AB and TWA 1/28 systems are physical binaries. 
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3.2. Optical Spectroscopy 

3.2.1. Spectral Morphology 



We classified TWA 30B by co mparison to M3-M5.75 rj Cha template s (denoted by RECX 



Lawson et al 



2001 



#) from iLuhman &: Steeghsl (|2004h , rj Cha has a similar age (~7 Myr; iMamaiek et al.l Il999l ; 



) to TWA; therefore, the templates should provide a better continuum fit with 
comparable surface gravities than would comparison to field templates, which have higher surface 
gravities. The TWA 30B spectra were boxcar smoothed down to comparable resolution of the 
templates (R ~1000) and then both data and template were normalized at 7500 A. The best match 
to the template spectrum over the wavelength range 6000-8000 A was determined by eye. This 
wavelength ra nge covers many broadba nd features typically used to classify M dwarfs such as TiO, 



VO, and CrH (jKirkpatrick et al 



1991 



)• 



The best fit to both epochs was to the M4 template RECX 5. We show the two unsmoothed 
epochs of MagE data along with the RECX 5 template in Figure [2j The continuum of TWA 30B 
shows little to no variation between epochs and matches well to the M4 template without having 
to account for variable reddening as is the case with TWA 30A (Looper et al. 2010). We therefore 
adopt an optical spectral type of M4 ± 0.5 for TWA 30B. 



3.2.2. Emission and Absorption Lines 

The continuum of TWA 30B varies little between the two epochs of MagE data. However, the 
forbidden emission lines (FELs): [S II], [O I], [O II], [O III], [N II], and [C I]; Na ID emission; and 
Li I and K I absorption are variable in both line strength and profile shape (Table U] and Figures [2] 
&[3|), while the Balmer emission lines are roughly constant (AEW < 15%). 

Like its comp anion, TWA 30B has FELs of [S II], [O I], and [N II] - signs of outflow and 



accretion in cTTS (jMuzerolle et al.l l2003). Notably, these lines are an order of magnitude stronger 
in the spectrum of TWA 30B than in TWA 30A. Particularly [N II] A6547, 6583, which are just on 
the threshold of detection in TWA 30A (EWs = -0.1 to -0.4 and -0.3 to -1.0 A, respectively) 
are enhanced to well-detectable levels in TWA 30B (EWs = -3.0 to -7.2 and -10.8 to -25.1 A, 
respectively). The strength of [O I] A6300 is increased from EWs of —3.1 to —15 A in TWA 30A 
to EWs of -85 to -226 A in TWA 30B, while the strength of [O III] A5007 is twice as strong 
in the spectrum of the latter source. This is a high excitation line that requires extremely high 
temperatures and low densities. 

All of the emission lines, with the exception of Ca II H+K, are stronger in the 2009 March 4 
spectrum of TWA 30B. The FELs increase by a factor of ~2-3 with the exception of [O III] A4958, 
5007 emission lines, which are fairly stable (AEWs = 0-30%; Table HJ). The gravity sensitive K I 
lines at 7665, 7699 A go into emission in the spectrum of 2009 March 4 while they are in absorption 
in the spectrum of 2009 May 30. This is also analogous to the changes seen in this line in TWA 
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30A. The Na ID A5890 feature in the 2009 May 30 data has a P Cygni profile, indicating outflowing 
winds along the line of sight. The Ha and H/3 lines are nearly constant in both line EWs (—7 A) 
and symmetric profiles in the two epochs. The full-width at 10% peak intensity of Ha is 220 km s _1 
and 200 km s _1 for the 2009 March 4 and 2009 May 30 data, respectively. The wid th of this line 
meets the >200 km s _1 threshold for accretion proscribed bv lWhite Basril (|2003l ). However, it 
is important to note that we expect this line to be b oth weakened and sy mmetrical if the disk is 
aligned edge-on, as is the case for MHO-5 in Taurus (jMohanty et al.1 12005). 



Another gravity sensitive line, Li I A6707, which is also used as a mass/age indicator (jBildsten et al 

19971). has a comparable stre ngth (EW = 0.46 ± 0.13 A) to other TWA M members (0.40 to 0.65 



Zuckerman & Sond 12004 ) in the 2009 May 30 data but is not detected in the 2009 March 4 



A, 

data. Together with the stronger emission line strengths in the spectrum of 2009 March 4, this 
indicates that the photospheric lines of Li I and K I are filled in (veiled) by enhanced accretion 
during this date. In Figure O we have plotted the radial velocity profiles of many of these lines 
after subtracting the radial velocity of TWA 30A (+12 km s _1 ). The centers of the photospheric 
lines K I A7699, Ha and H/3 match this stellar rest velocity (A v < 5 km s _1 ). All other lines 
(FELs, emission, and absorption) show only small deviations from the stellar rest frame (A v < 30 
km s _1 ; Table U]), suggesting that associated outflows are observed in the plane of the sky, so that 
the (assumed perpendicular) circumstellar disk is viewed edge-on. 

TWA 30B shows no He I emission at 4471, 5876, 6678, and 7065 A. He I is present in the 
accreting late-M members TWA 2 7AB and TWA 28 a nd weakly in emission in the accreting late-K 



TWA 1 and the early-M TWA 3 (IHerczeg et al J 120091 ). On the other hand, the spectrum of TWA 



30B exhibits [C I] A8727, 9824, and 9850 in emis sion in both ep o chs, w hich is not present in the 
spectrum of TWA 30A or any known PMS stars. lErcolano et al.l (|2009l ) have suggested the use of 
this line in combination with the fine structure MIR lines of [C I] to constrain the size of gaseous 
protoplanetary disks. Their model for a T e ff = 4000 K, log g = 3.5, R = 2.5 R star predicted 
EWs of 2-5 mA for the NIR [C I] lines. We have measured EWs of -0.4 to -3.4 A: a factor of 
1000 times greater in strength. 



3.3. NIR Spectroscopy 

3.3.1. Spectral Morphology 

In the left panel of Figure HI we show seven NIR spectra of the 18 epochs, which span the 
full range of colors and dates observed. Since we did not have a set of young M dwarfs observed 
at similar wavelengths and resolution as our data, we classified the TWA 30B NIR spectra by 
comparison to field M dwarfs taken with an identical set-up (Kirkpatrick et al. 2010): Gl 752A 
(M3V), Gl 213 (M4 V), Gl 51 (M5V), LHS 1375 (M6 V), and vB8 (M7 V). In the right panel of 
Figure HI we show the first epoch of our TWA 30B NIR data from 2009 May 14 UT in comparison to 
the M4 standard Gl 213. The strengths of the TiO and H2O bands as well as the overall morphology 
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from 0.7-1.3 /xm match well between the standard and data for all 18 epochs. This region is similar 
to the wavelength coverage (0.9-1.4 /im) proscribed by Kirkpatrick et al. (2010) for NIR spectral 
classification. At longer wavelengths (H and K), the spectra of TWA 30B showed variable levels 
of NIR excess. We therefore adopt a NIR spectral type of M4 ± 0.5. 

We measured the spectrophotometry of these data using the 2MASS JHK S filter profiles 
and report these values in Table [5] and plot them in JHK S color-color space (along with TWA 
30A) in Figure [5j These epochs have J — K s measurements of 1.1-2.2 mag, much redder than 
the 0.86 ± 0.16 m ag average of 2126 field M4 dwarfs in the SDSS DR5 sample with 2MASS colors 
((West et al-lbood ). The original 2MASS J-K s color is 1.63 ± 0.06 mag. TWA 30B therefore shows 
a variable amount of NIR excess, with a sharp increase from 2009 December 28 UT (J — K s = 1.3 
mag) to the 2010 January 2-3 UT epochs (J — K s = 1.9 to 2.2 mag) before settling down to a 
more modest NIR excess in the 2010 January 24-28 UT data (J — K s = 1.1 to 1.2 mag). Finally 
in the 2010 May data, the NIR excess increases again (J — K s = 1.6 to 2.0 mag). Gr aphically, we 



see t hat this NIR color variation in TWA 30B is nearly parallel to the cTTS locus (IMeyer et al 



1997^ in Figure [5l w hile the NIR variability in TWA 30A is parallel to the reddening vector from 
Cardelli et al.l (|1989i ). The slope differences in color tracks suggests that the excesses in the these 
two sources arise from different origins or viewing geometries. 



3.3.2. Emission and Absorption Lines 

Most of the spectroscopic NIR data we obtained for TWA 30B is too low resolution (R ~150) to 
examine the emission and absorption line features. We note, however, that even at these resolutions 
strong He I emission appears at 1.083 /tin and blended [C I] emission at 0.982 and 0.985 /mi. We 
have verified the wavelengths of both of these features using our three epochs of moderate resolution 
(R ~1200) NIR SXD data. In Figure SI we show the SXD data obtained on 2010 Jan 3, covering 
the emission lines of He I and [C I] in the J-band. 



4. Discussion 

4.1. The Nature of NIR Excess Emission in TWA 30B 

As noted, the shape of the 0.4-1.4 /mi spectrum of TWA 30B is relatively steady and consistent 
with an M4 spectral classification while the 1.4-2.5 /mi spectrum, in contrast, is highly variable. 
This dichotomy suggests that the long wavelength flux originates from a distinct and variable 
component. To isolate this NIR excess, we subtracted the M4 template Gl 213 from each epoch of 
data and fit a black body (BB) curve to the residuals. BB fluxes were calculated for temperatures 
spanning 500-1500 K (in 10 K increments) and initially scaled to a source radius of 1 R© and 
a distance of 42 pc. For each spectrum, we renormalized the BB fluxes to match the residual 
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flux at 2.2 /jm. To incorporate uncertainties in these spectral observations, we repeated the fits 
by adding a random Gaussian noise component pixel-by-pixel based on the noise spectrum. This 
procedure was done 100 times with the set of BB curves fit to each iteration. We then calculated 
the reduced x 2 values between the set of BB curves and the residuals over the wavelength ranges 
of 1.45-1.8 and 2.0-2.35 /im, excluding the region between these two bands affected by telluric 
absorption. The best-fitting BB curve for this fit and the associated normalization were adopted as 
the temperature and visible area (in units R© 2 ), respectively, of the emitting region. The adopted 
parameters (temperature and area) and uncertainties for each epoch are the averages and standard 
deviations for these simulated datasets. These values are listed in Table [5} an example of a BB 
best-fit is shown in the right panel of Figure HI 

To interpret our results, we make a few assumptions about the origins of the infrared excess. 
Given the good fit to a single temperature blackbody, we assume that it is reprocessed stellar 
radiation, emitted by some small (approximately single temperature) region of the dislJE The BB 
temperature observed will then depend on both the radial distance from the star to the emitting 
region and the inclination of the region to the stellar photons. Lacking a priori knowledge of the 
inclination, we first consider the two limiting cases: (1) an annulus in a thin flat disk (i.e., the 
emitting region is confined to the plane of the disk), and (2) a vertical disk rim or warp (i.e., the 
emitting region is perpendicular to the plane of the disk, therefore normal to the incident stellar 
radiation). For a star of radius R* and effective temperature T#, the temperature of a disk 
region at a radial distance Rd for case (1) is given by: 

rp4 1 * 



which holds in the limit Rd ^> R* (see Equation 
For case (2) the analogous expression is: 



2 fR£V 

3 \Rd) 



[1] 



[A22] in lAdams et al.lll98 



T 4 = T 4 

Both equations contain the usual dilution factor 
having an additional projection factor that goes 



(I) 2 ' 2 > 

of (R*/Rdf. The thin flat disk case differs in 
to R*/Rd for large Rd (see footnote [4]). For 



3 As such, we neglect the possibility that it is simply starlight scattered by the disk, since to zeroeth order this 
would not yield an excess spectral energy distribution mimicking a single temperature BB. 

11988), retaining only the reprocessing term - reads: T^ 4 = 
(R,/Rd)- Using a Taylor expansion in the limit of u < 1 gives 



4 The exact equation 



[A22] 

2M/21 



Adams et al 



(T* /7r)[arcsin(it) — u(l — u ) ' ], where u 
T d A = (T» 4 /7r)(2u 3 /3), which is our equation [1] above. Comparing the value of [arcsin(u) — u(l — u 2 ) 1 ^ 2 ] to that of 
its Taylor expansion approximation (2u /3), results in a difference of <5% for u < 0.4; i.e., the use of the approxi- 
mation in equation [1] is valid, with very small errors, for any Rd > 2.5R*. 
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inclinations intermediate between a flat disk and a vertical rim/warp, the disk temperature will 
vary approximately as Td oc R^ 9 , with q ranging between the limiting values of 3/4 (eqn. [1]) and 
1/2 (eqn. [2]). Conversely, for a given BB temperature Td, equations [1] and [2] provide the lower 
and upper limits, respectively, on the radius Rd from which the emission can arise. 

Knowing the BB disk temperatures from our fits to the observed excess, we invert equations 
[1] and [2] to estimate the corresponding Rd limits. We assume T* = 3200K and R* = 1 Rq, 
appropriate for a youn g M4 star, using the s pectral type to effective temper ature conversion scal e 
for young objects from iLuhman et al.l (|2003l ) and evolutionary models from iBaraffe et al.l ([1.998). 
These adopted values correspond to a mass of ~0.15 Mq at an age of ~10 Myr. The inferred Rd 
are listed in Table We see that, for a thin flat disk, Rd is of order 3-5 i?*, while for a vertical rim 
or warp, the values are a factor of ~5 higheJl. We note that the small Rd values calculated for the 
thin flat disk case are consistent with the inner edge of a magnetosphericaUy-truncated accretion 
disk, where disk structure is expected (see below). 

In both cases, a lower Td is associated with a larger Rd, which follows from our assumption 
that the radiation is reprocessed starlight from a disk region with a fixed inclination (either from 
a flat disk or a vertical rim/warp). However, the Td variations may equally well arise from re- 
gions with varying inclination at roughly constant Rd, i.e., by changing q instead of Rd in the 
general reprocessing approximation Td oc R d 9 ■ For example, a rim or warp can f o rm at the inner 



edge of the disk, due to frontal heating by stellar irradiation (jNatta et al 



2001 



D'Alessio et al 



2005 ) and/or could be formed from the interaction between the stellar magnet o-sphere and the dis k 



(jTerquem h Papaloizou 



2000 



Bouvier et al.ll2003l : also invoked for TWA 30A: lLooper et alJl201tJ ). 



So suppose we assume that the emission arises from a fixed radius at the disk inner edge, with the 
lowest observed Td arising from a flat disk geometry and higher values arising from a burgeoning 
warp, i.e., the emission region at larger inclinations. Table [5] shows that our smallest Td is consis- 
tent with flat disk emission from Rd ~ 5 R*. For Rd fixed at this value, the highest Td would then 
require an inclination corresponding to q ~ 0.6, and the remaining Td would require 0.6 < q < 0.75: 
in all cases intermediate between the flat and vertical values of 0.75 and 0.5, respectively. 

So far we have discussed the geometrical constraints - location and inclination - on the emitting 
region implied by our inferred BB temperatures, assuming the excess flux is reprocessed stellar 
radiation from the disk. In the next section (4.3), we present several lines of evidence suggesting 
that the disk around TWA30B is edge-on to our line of sight. Therefore, the temporal variability in 
the NIR excess may be most straightforwardly interpreted as spatial variations in the disk structure 
rotating into our line of sight and out of it (occulted by the star and foreground parts of the disk). 
If so, relatively smooth azimuthal variations in disk structure should cause variability on timescales 
comparable to the orbital period at the emitting radius. For the stellar mass and radius assumed, 
the Rd ~ 3-5 i?* derived for the thin flat disk case (and corresponding to the expected location 



5 Note that the Rd inferred for the thin fiat disk case are all large enough - Rd > 2.5R* - to be a posteriori 
consistent with the use of the approximate equation [1], as explained in footnote [4]. 
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of the disk inner edge) imply Keplerian orbital periods of ~2-4 days (Table [5]), consistent with 
the observed variations in the excess over timescales of a day. The Rd ~ 15-25 i?* derived for 
a purely vertical rim/warp, on the other hand, imply orbital periods of ~15-40 days (Table [5]), 
much longer than the observed timescales. This suggests that the inner edge of the disk, where 
disk structure may naturally arise as noted above, may indeed be the source of the NIR excess. We 
note that if an edge-on disk is significantly occulting the star (as argued in Section 4.3), then we 
are unlikely to be directly observing the inner regions of the disk where the NIR excess appears to 
arise (since the stellar occultation is more plausibly caused by the flared outer regions of the disk, 
which would occlude the inner disk as well). In this case, it is likely that the inner disk emission is 
being scattered into our line of sight by surrounding disk regions (see also below). 

Lastly, we consider trends in the area of the emitting region inferred from our BB fits. Fig. 
6 shows the derived source areas plotted against the BB temperatures for all the observed SpeX 
prism spectra. It is clear that the emitting area decreases with increasing temperature. To explain 
this trend, we use the following salient points: the areas we derive do not correspond to the total 
area of the emitting region, but to the fraction of that area from which radiation actually reaches 
us; the variability timescale of the excess is consistent with emission from the disk inner edge; the 
derived Td are consistent with reprocessed emission from a region at fixed Rd corresponding to this 
edge, but with variable inclination; and the observed emission is likely not reaching us directly, 
but via scattering off surrounding disk surfaces. Now, for a fixed Rd, the associated Td increases 
with increasing inclination (i.e., as the source region becomes more vertical). Thus the inverse 
trend in area with Td can be interpreted as a decrease in the apparent source area with increasing 
inclination, if indeed Rd is roughly constant. A more vertical region will emit a larger fraction 
of its radiation towards the central star, leaving less to be scattered by surrounding disk surfaces 
into our line of sight. This would automatically yield a reduction in the apparent source area with 
increasing inclination (and hence with increasing Td for fixed Rd), as observed. 

The scenario postulated above - reprocessed starlight emitted by disk structure close to the 
disk inner edge, and variably scattered towards us - is clearly a simplistic toy model: while it is 
consistent with the temperatures, apparent source areas and variability timescales of the observed 
NIR excess, it is only loosely constrained so far by these parameters. Further observations at higher 
temporal sampling, and extended to longer wavelengths are required to establish its veracity. 

4.2. Comparison of TWA 30 A and 30B 

Both components of this very wide (p ~3400 AU) binary share many physical properties, such 
as a near-equivalent spectral type (M5 ± 1 for TWA 30A and M4 ± 0.5 for TWA 30B) and strong 
FELs of [O I], [O II], [N II], and [S II]. These FELs are spectroscopic signatures of an outflow. 
In both cases the small velocity shifts (<30 km s _1 ) of these lines from the rest velocity of the 
star suggest that the outflow lies perpendicular to the line-of-sight and, by extension, that the disk 
orientation lies close to edge-on. The narrow P Cygni feature in the Na ID lines of both stars 
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also points toward the outflow being viewed perpendicular to the line of sight. The presence of 
[O III] in the spectrum of cTTS is rare an d has only been prev i ously reported in the literature in 
two high accretion ra te cTTS - DG Tau (jHartigan et all Il995l ; feaga et al.1 Il996l ) and XZ Tau B 



([White k Ghedl200ll ) - and i n TWA 30A (Lqoper et al. 2 010). This high-ionization line indicates 



Schwartz 


1983; 


Bohm 


1995) 



collisional shock fronts (e.g., 
winds from massive stars, where the requisite extremely high temperatures and l ow densities can 



be fo und. Given that the TWA 30 AB system resides far from any molecular cloud ([Tachihara et al 



2009), Looper et al. (2010) posit that this line may be formed by newer outflow material overtaking 
and collisionally shocking older jet fronts. 

With strong spectral signatures of outflows, it is at first surprising that the typical signatures of 
accretion in cTTS - strong Ha (usually defined as > — 10 A) and He I emission at 4471, 5876, 6678, 
and 7065 A - are weak or ab sent in TWA 30A and TWA 30B. The EWs of Ha for both components 
is —7 A. lLooper et al.1 (2010j) posit that Ha emission, which in addition to chromospheric origins is 
formed in the shock region where infalling accreting gas impacts the stellar surface, is obscured by 
optically thick gas from a nearly edge-on disk. He I is also formed in the shock region, producing 
the narrow component, and in the accretion flow, producing the broad component. TWA 30A was 
weakly detected in the RASS Faint Source Catalog as having X-ray emission while no such detection 
was made for TWA 30B. An edge-on configuration could likewise obscure the X-ray emitting region 
of the star. 

Many of the emission lines, particularly FELs, are common to spectra of both stars, but a few 
differences exist. Mg I] A4571 emission appears in the spectrum of TWA 30A but not in TWA 30B. 
The detection of this line is similar in rarity to [O III] emission in cTTS and is seen in o nly two 
other PMS stars - XZ Tau B, which also has [O III] emission, and the EXor VY Tauri (M0; iHerbie 
1990). This line indicates the presence of gas in the circumstellar disk. The detection of Mg I] in 
higher mass cTT S remains prohibitive ly difficult due to the strong continuum from the star around 
this wavelength (jErcolano et al.l l2008). 



[C I] emission is present in both our optical and low-resolution NIR spectra in T WA 30B but is 
unseen in TWA 30 A. This line has been d etected in the spectra of pla n etary nebulae (Jewitt et all 



1983 



Liu et al.1 1 19951 ) , mo lecular clouds (IHippelein fc Muenchl Il978l: iMunch fc HippeleirJ Il982l ) 



R Coronae Borealis stars (jPandev et al.1 120041 ). and cometary tails (jOliversen et al.1 120021 ). The 
presence of [C I] in the spectrum of a PMS star has not been previously reported in the literature. 
The ratio of this line with the fine structure MIR lines of [C I] offers a potentially powerful tool to 
constrain the size of the gaseous component of a protoplanetary disk (jErcolano et al.ll2009i ). As of 
yet, no mid-infrared (MIR) observations of TWA 30B have been made. 

TWA 30A and TWA 30B present similar lines of evidence for a near edge-on orientation of 
an accreting circumstellar disk. However, it is apparent that the disk of TWA 30B is viewed fully 
edge-on and is occulting the star. This is highlighted from the fact that TWA 30B appears ~5 
mag underluminous in optical and NIR bands compared to TWA 30A despite being one spectral 
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type earlier (Table [3|). TWA 30A has temporal c hanges in both emis sion line strengths profiles 
and variable reddening (Ay = 1.5-9.0 mag), leading lLooper et al.l (j2010n to suggest that the stellar 
axis may be inclined to the plane of the disk. While we do not have the same coverage in the 
optical data of TWA 30B as obtained for TWA 30A, the two optical spectra of TWA 30B and the 
J-band portion of the 18 epochs of low-resolution NIR data show little variation in spectral slope or 
evidence for reddening. Instead, the NIR excess present in the H- and -fT-band portions of the NIR 
spectral data show large variations between epochs on timescales as short as one day. In Section 
4.1, we proposed that these changes may arise from disk structure, such as a rim or a warp, at the 
inner disk edge that emits reprocessed starlight. 

Assuming a distance of 42 ± 2 pc, TWA 30A and TWA 30B would be the 2nd and 3rd closest 
accretors to t he Sun. The closest accreting star is anothe r TWA system, the M4 PMS star TWA 
3 (Hen 3-600; lHenizdll976l ). with a distance of 34 ± 4 pc (IMama j ekl 120051 ) . Distances to the other 
nearby accreting TWA members are: TW Hya (53 ± 3 pc), TWA 28 (55 ± 1.5 pc), and TWA 27AB 
(53 ± 1 pc), from Mamajek, in prep. Both stars in the TWA 30AB system power the strongest 
jets of any TWA systems based on the strengths of their FELs. Given the proximity of the system, 
they are excellent candidates for follow-up imaging to resolve their outflows. 



4.3. Are TWA 30AB Members of the f3 Pictoris Group? 



While TWA 30A and TWA 30B have kinematics and spectroscopic indicators pointing towards 
membership in the TWA, there exists another notable young group of s tars within tens of pc of 
the Sun - the f3 Pic moving group (~12 Myr; IZuckerman et al.1 l2001al ). There has been some 
recent discussion on whether some previously clai med TWA stars c ould actually be j3 Pic group 
members. Most notably, the low-mass star TWA 22 (jSong et al.ll2003) was claimed to be the closest 
TWA member. It was subsequently rejected as a TWA member by iMamajekl (120051 ) based on its 
poor fit in the conver g ent p oint solutio n to the group . Thi s claim of non-membership was then 
refuted by ISong et al.l (|2006l ) . Recently, iTeixeira et al.l (|200 9) showed that TWA 22 is much more 
likely to be a f3 Pic member than a TWA member based on a precise space velocity and kinematic 
traceback, enabled by a precise trigonometric parallax. The tale of the membership of "TWA" 22 
reiterates that the space motions of nearby young groups are fairly similar at the km s _1 level, and 
membership may not be obvious. 



Using the space motion for the f3 Pic group from lZuckerman &: Songl (J200J), we compared the 
proper motions and radial velocities for TWA 30A and TWA 30B to that expected for an ideal f3 
Pic member. The peculiar (perpendicular tangential) velocities for TWA 30A and TWA 30B would 
be 3.0 ± 0.3 km s _1 and 1.6 ± 1.1 km s , respectively, and the predicted radial velocities would 
be 6.9 km s _1 each if (3 Pic members. In comparison, if TWA 30A and TWA 30B are members of 
the TWA, the peculiar velocities are only 1.2 ± 0.3 km s _1 and 0.1 ± 0.9 km s _1 , respectively. The 
measured radial velocities of TWA 30A and TWA 30B are 12.3 ± 1.5 km s" 1 and 12 ± 3 km s -1 , 
respectively; well in-line with the predicted radial velocity of 11.0 db 0.9 km s _1 for each if they are 
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true TWA members. Given the kinematic data in hand, we reject the hypothesis that the TWA 30 
pair could be f3 Pic members rather than TWA members. 



5. Conclusion 



We report the discovery and characterization of the M4 PMS star TWA 30B, which has a 
projected separation of ~3400 AU from the M5 PMS star TWA 30A. This binary is a member of 
the TWA and has stronger spectroscopic signatures of outflows than any other TWA members, with 
FELs of [O I], [O II], [N II], and [S II]. Both members display strong [O III] emission, which is rare in 
cTTS and possibly indicates collisional shock fronts from the jets. All FELs have negligible radial 
velocity shifts from the stellar rest frame for both TWA 30A and 30B, suggesting that the outflows 
from both components are oriented in the plane of the sky. TWA 30B is underluminous by ~5 mag 
in both optical and NIR photometric bands in comparison to TWA 30A, suggesting that an edge-on 
disk heavily obscures the stellar photosphere. TWA 30B has variable NIR excess well-fit by black 
body fluxes, with temperatures ranging from 630 to 880 K and emitting areas that are inversely 
proportional to the temperature. The detection of [C I] emission at 8727, 9824, and 9850 A mark 
the first detection of this line in a PMS star. If paired with a future measurement of the [C I] 
fine-structure lines in the MIR, these measurements could afford an independent measurement of 
the gaseous component of the circumstellar disk. The TWA 30AB system represents an invaluable 
case study in low-mass stars as the 2nd and 3rd closest accretors to the Sun (42 ± 2 pc), marked by 
a wide array of forbidden emission lines and temporal variability in NIR excess arising from each 
star's circumstellar disk. 
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Table 1. Astrometry for TWA 30B 



RA (deg) 


Dec (deg) 


a RA (mas) 


a Dec (mas) 


Epoch (yr) 


Survey 


Ref 


173.076795 


-30.308531 


200 


200 


1968.7 


USNO-B1.0 


1 


173.076125 


-30.308692 


207 


200 


1994.1 


USNO-B1.0 


1 


173.075928 


-30.308792 


80 


60 


1999.23 


2MASS 


2 


173.075736 


-30.308848 


20 


20 


2008.95 


ISPI 


3 


173.075692 


-30.308861 


20 


20 


2009.19 


ISPI 


3 


173.075672 


-30.308865 


20 


20 


2009.92 


ISPI 


3 


173.075711 


-30.308872 


20 


20 


2010.08 


ISPI 


3 



Note. - R eferen ces: (1) USNO-B1.0 (Monet et al.1 120031 ). (2) 2MASS (jCutri et al.1 12003 : 
Skrutskie et al.ll2006l ) and (3) this paper. 
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Table 2. Spectroscopic Observation Log for TWA 30A and 30B 



W LJJ L 


Tel /Tnst 




A/AA 


TIT Datp a 


N x t (i\ b 


Z c 


Calibrator 


v_-'U1ilii Li una 


TWA 30B 


Magellan/MagE 


0.30-1.05 


4100 


090304 


1 x 1800 


1.0 


Hiltner 600 (WD) 


Clear, 0'.'7 seeing 


TWA 30B 


Magellan/MagE 


0.30-1.05 


4100 


090530 


1 x 2400 


1.1 


GD 108 (WD) 


Thin Clouds, 0'.'9 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


090514 


8 x 150 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'5 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


090628 


4 x 120 


2.3 


HD 98949 


(AO V) 


Clear, 0'.'5 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


090629 


4 x 180 


2.7 


HD 98949 


(AO V) 


Clear, 0'.'6 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


091228 


10 x 150 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'6 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100102 


10 x 150 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'8 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100103 


10 x 150 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'6 seeing 


TWA 30B 


IRTF/SpeX 


0.80-2.40 


1200 


100103 


24 x 180 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'6 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100124 


12 x 150 


1.6 


HD 98949 


(AO V) 


Light Cirrus, 0''6 seeing 


TWA 30B 


IRTF/SpeX 


0.80-2.40 


1200 


100124 


26 x 180 


1.6 


HD 98949 


(AO V) 


Light Cirrus, 0'.'6 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100125 


16 x 150 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'7 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100126 


20 x 150 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'9 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100127 


10 x 150 


1.9 


HD 98949 


(AO V) 


Clear, 0'.'9 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100128 


18 x 150 


1.7 


HD 98949 


(AO V) 


Clear, 0'.'9 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100520 


10 x 120 


1.6 


HD 98949 


(AO V) 


Clear, 1'.'5 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100521 


12 x 180 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'7 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100522 


6 x 180 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'5 seeing 


TWA 30B 


IRTF/SpeX 


0.80-2.40 


1200 


100522 


24 x 180 


1.7 


HD 98949 


(AO V) 


Clear, 0'.'5 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100523 


10 x 180 


1.8 


HD 98949 


(AO V) 


Clear, 1'.'5 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100525 


10 x 180 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'6 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100526 


10 x 180 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'7 seeing 


TWA 30B 


IRTF/SpeX 


0.70-2.50 


150 


100527 


10 x 180 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'7 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


091229 


10 x 120 


1.6 


HD 98949 


(AO V) 


Clear, 0''4 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


091231 


16 x 120 


1.7 


HD 98949 


(AO V) 


Clear, 0'.'3 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


100102 


6 x 120 


1.5 


HD 98949 


(AO V) 


Clear, 0'.'8 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


100124 


6 x 120 


1.8 


HD 98949 


(AO V) 


Clear, 0'.'6 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


100125 


22 x 120 


1.7 


HD 98949 


(AO V) 


Clear, 0'.'7 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


100126 


12 x 120 


1.8 


HD 98949 


(AO V) 


Clear, 0'.'8 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


100520 


10 x 120 


1.6 


HD 98949 


(AO V) 


Clear, 1'.'5 seeing 


TWA 30A 


IRTF/SpeX 


0.80-2.40 


1200 


100525 


10 x 120 


1.6 


HD 98949 


(AO V) 


Clear, 0'.'9 seeing 



a Epochs are denoted as YYMMDD in all tables. 
b Number of integrations times the integration time. 
c Airmass. 
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Table 3. Properties of TWA 30 AB 



Parameter TWA 30A TWA 30B Adopted Mean Ref 



2MASS Designation 3 - 


J11321831 -3019518 


Jl 1321822— 3018316 




1 


a°(J2000) a 


173.076322 


173.075928 




1 


<5°(J2000) a 


-30.331059 


-30.308792 




1 


/x Q (masyr L ) 


-89.6 ± 1.3 


-83 ± 9 




2,3 


/ — 1 \ 

(masyr L ) 


-25.8 ± 1.4 


-30 ± 9 




2,3 


Distance (pc) 


42 ± 2 


44 ± 5 


42 ± 2 


3,4 


\ rad (km s x ) 


12.3 ± 1.5 


12 ± 3 


12 ± 2 


3,4 


A /TV T \ 

Age (Myrj 




~8 


~8 


3,4 


Optical SpT 


M5 ± 1 


M4 ± 0.5 




3,4 


±> (mag) 


15.6 


20.4 




5 


it (mag) 


12.9 


17.4 




5 


1 (mag) 


11.30 ± 0.03 


16.68 ± 0.08 




6 


DENIS J (mag) c 


9.78 ± 0.06 


14.00 ± 0.09 




6 


DENIS K (mag) c 


8.88 ± 0.07 


11.08 ± 0.08 




6 


DENIS J-K (mag) c 


0.90 ± 0.09 


2.92 ± 0.12 




6 


2MASS J (mag) c 


9.64 ± 0.02 


15.35 ± 0.05 




1 


2MASS H (mag) c 


9.03 ± 0.02 


14.53 ± 0.05 




1 


2MASS K s (mag) c 


8.77 ± 0.02 


13.72 ± 0.04 




1 


2MASS J — K s (mag) c 


0.87 ± 0.03 


1.63 ± 0.06 




1 


Angular Separation 13 (") 






80.2 


3 


Projected Separation 13 (AU) 






3370 ± 160 


3 


Est. Actual Separation 13 (AU) 






4250 ± 200 


3 


Li EW (A) 


0.6 ± 0.1 


0.5 ± 0.1 d 




3,4 


Ha EW (A) 


-6.8 ± 1.2 


-7.4 ± 0.02 e 




3,4 



a From the 2MASS Point Source Catalog at epoch 24 March 1999 UT. 
b Estimated based on a convergent point analysis (see Section 3.1). 

C TWA 30A and TWA 30B have time variable photometry (Looper et al. 2010 and this paper). 

d Derived from the equivalent width measured from the MagE data of epoch 2009 May 30 reported 
in Table H We did not detect Li I in the MagE data of epoch 2009 March 4. 

e Derived from the average equivalent widths measured from the two MagE spectra reported in Table 
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II 



Note. — References: (1) 2M ASS (jSkrutskie et alJl2006h. (2) UC AC3 (jFinch et al.1 201Gh. (3) this 



pape r, (4) lLooper etaD (|201ol l. (5) USNO-A2.0 (jMonet et al1ll998h . and (6) DENIS (|Epchtein et al 
19971 ). 
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Table 4. Equivalent Widths for Selected Line Features of TWA 30B 



UT 


\ab (A) 


Kbs (A) 


Ion 




Flux a 


EW ± Iff (A) b 


V (km s _1 ) c 


090304 


3934 777 


3935 06 


Ca II K 






—24.5 


22 


090^30 


3934.777 


3934.93 


Ca II K 






— 104 


12 


090304 




3969 90 




- He d 




—3.3 


23 


uyuoou 




oyuy. f o 


P a IT FT 4 


H^ d 
-lit 




25 1 




090304 


4069 749 


4070 029 


fS III 




66.2 


—61 ± 10 


21 


uyuoou 


4f)fiQ 74Q 




ra Til 
[5 11] 






— 19 2 




090304 


4077 500 

rtU I 1 . UUU 


4077.70 


[S Til 

[O 11] 






— 14.5 


15 


090^30 

uyuoou 


4077 ^oo 

'iU 1 f . UUU 


4077.73 


fS TT1 

[3 11] 






— 1.9 


17 


090304 


4862 683 


4862 935 


H/S 




11.4 


—8 ± 2 


16 


090^30 

uyuoou 


4869 683 
iou^ . UOO 


4869 90 
iouz . yu 


up 






— 7 ± 2 


13 


090304 


4960 son 


4960 765 


ro mi 

[VJ 111] 




24.9 


— 19 ± 1 


28 


090^30 

uyuoou 


4960 300 
iyuu . ouu 


4960 70S 
rtyuu. 1 uo 


fO TTT1 

[VJ 111] 




24.8 


— 19 i 3 


24 


090304 


5008 240 


5008 465 


TO TTT1 

[VJ in] 




103 


—94 ± 2 


13 




5008 240 


5008 755 

JUUO. 1 OO 


[O TTT1 

[VJ 111] 




81.2 


— 72 ± 3 


31 


090304 


5200.621 


5200 438 


\Fc TT1 
[re iij 




11.9 


—7 ± 2 


_11 


090^30 


^900 691 

OirUU . J- 




[Fo 111 
[fe 11] 






<r —4 




090304 


5578 887 


5579 01 3 


[O 11 

[VJ 1] 




11.0 


—6.5 ± 0.4 


7 


090^30 

uyuoou 


^78 887 


^79 000 
oo / y.uuu 


[O Tl 
[U 1J 




4.6 


— 1.4 ± 0.6 


Q 


090304 


5756 240 


5756 678 


[N Til 




10.6 


—4.2 ± 0.7 


23 


090^30 

uyuoou 


^7^6 940 




\N TT1 

[IN 11J 










090304 


5891 583 


5892 409 


Na I D 




21.3 


— 18.1 ± 0.8 


42 


090^30 

uyuoou 


^891 ^83 


^891 6^ 
ooy i .uo 


Na I D 






—3.8 


3 


090304 


5897 558 


5898 340 


Na I D 




18.2 


— 15.1 ± 0.8 


40 


090530 


5897 558 


5897 67 

OOC 1 .U 1 


Na I D 






—2.6 


g 


090304 

Uf UOUrt 


6302 050 


6302 283 


[O Tl 

[VJ 1] 




233 




H 


090^30 

uyuoou 


6309 0^0 
uou^ . uou 


6309 1 98 
uouz . -L yo 


[O Tl 

[VJ 1J 




92.0 


_Qzl 7 4- n 6 

— O^l. 1 ^ U.U 


7 


090304 


6365 540 


6365 789 


ro ti 

[VJ 1] 




55.8 


_r;n i _|_ n o 

OU. J. _1_ U . >J 


12 


VJ fU'JOU 


6365 540 

UOUO . <JrtU 


6365 696 

UOUO . VJ CVJ 


ro ti 

[W 1J 




26.5 


_20 9 -1- 1 

ZjXJ . tJ _1_ U . -L 


7 


000304 


6549 850 


6550 1 07 


r\ i ii 

[IN 11] 




13.0 


— 7.2 ± 0.2 


12 


uyuoou 


DD4y .ODU 


OOOU.U f 1 


T1\T TTl 

[JN 11J 




o.y 


on I rt -i 
— O.U ± U.l 


1U 


090304 


6564.610 


6564.779 


Ha 




14.4 


-8.1 ± 0.2 


8 


090530 


6564.610 


6564.738 


Ho 




13.3 


-7.2 ± 0.6 


6 


090304 


6585.280 


6585.526 


[N II] 




31.5 


-25.1 ± 0.2 


11 


090530 


6585.280 


6585.417 


[N II] 




17.6 


-10.8 ± 0.1 


6 


090304 


6709.660 




Li I 






<0.2 




090530 


6709.660 


6709.831 


Li I 




3.0 


0.46 ± 0.13 


8 


090304 


6718.290 


6718.528 


[S II] 




10.2 


-5.8 ± 0.2 


11 


090530 


6718.290 


6718.404 


[S II] 




7.4 


-2.9 ± 0.2 


5 


090304 


6732.670 


6732.909 


[S II] 




16.4 


-10.4 ± 0.2 


11 


090530 


6732.670 


6732.891 


[S II] 




11.4 


-6.0 ± 0.2 


10 


090304 


7322.010 


7322.258 


[O II] 




18.8 


-10.7 ± 0.1 


10 


090530 


7322.010 


7322.316 


[O II] 




12.7 


-4.7 ± 0.1 


13 


090304 


7331.680 


7332.663 


[O II] 




15.4 


-7.8 ± 0.1 


40 


090530 


7331.680 


7332.495 


[O II] 




10.5 


-3.4 ± 0.1 


33 


090304 


7667.021 


7667.435 


K I 




14.5 


-6.1 ± 0.4 


16 
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Table 4 — Continued 



UT 


Mab (A) 


Kbs (A) 


Ion 


Flux a 


EW ± la (A) b 


V (km s- 1 ) 


090530 


7667.021 




K I 








090304 


7701.093 


7701.61 


K I 




-3.1 


20 


090530 


7701.093 


7701.335 


K I 


3.5 


0.46 ± 0.07 


9 


090304 


8185.505 


8185.537 


Na I 


4.1 


0.76 ± 0.12 


1 


090530 


8185.505 


8185.800 


Na I 


3.6 


0.66 ± 0.09 


11 


090304 


8500.36 




Ca II 




<0.2 




090530 


8500.36 




Ca II 




<0.2 




090304 


8544.44 


8545.08 


Ca II 




0.52 


22 


090530 


8544.44 


8545.32 


Ca II 




0.86 


31 


090304 


8664.52 


8664.73 


Ca II 




0.62 


7 


090530 


8664.52 


8665.07 


Ca II 




0.73 


19 


090304 


8729.520 


8729.531 


[CI] 


11.4 


-2.4 ± 0.1 





090530 


8729.520 


8729.865 


[C I] 


6.8 


-0.4 ± 0.1 


12 


090304 


9826.820 


9826.803 


[CI] 


5.6 


-1.3 ± 0.1 


-1 


090530 


9826.820 


9827.020 


[C I] 


4.9 


-0.6 ± 0.1 


6 


090304 


9852.960 


9853.173 


[CI] 


10.7 


-3.4 ± 0.1 


6 


090530 


9852.960 


9852.925 


[C I] 


6.4 


-1.4 ± 0.1 


-1 



a The integrated line fluxes are given in units of 10 — 16 erg cm -2 s~ 1 and should only be 
used to calculate relative line fluxes between features in the same spectrum as our data are 
not photometrically calibrated and hence do not account for slit losses or non-photometric 
conditions. 

b Some lines were measured by hand and have typical uncertainties of 0.5-3 A. 

c Reported as heliocentric velocities, i.e., not relative to the assumed stellar RV of 12 ± 3 
km s _1 . 

d This line is blended or contaminated. 
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Table 5. Derived NIR Magnitudes 51 , Colors 3, , and Disk Parameters of TWA 30 AB 



Object UT b J c J-H c H — K s c J - K a c BB (T d ) d Area 6 R<n f ' 9 Pl l ' h Rd2 f ' 9 P2 f ' h Notes 



TWA 30B 


990324 


15.35 


0.82 


0.81 


1.63 


















2MASS PSC 


TWA 30B 


090514 


15.4 


0.70 


0.69 


1.39 


670 


± 


20 


4.4 ± 0.7 


4.8 


3.2 


22.8 


33.0 


SpeX prism 


TWA 30B 


090628 


15.3 


0.72 


0.50 


1.22 


880 


± 


40 


0.45 ± 0.09 


3.3 


1.8 


13.2 


14.6 


SpcX prism 


TWA 30B 


090629 


15.3 


0.73 


0.56 


1.29 


860 


± 


40 


0.48 ± 0.10 


3.4 


1.9 


13.8 


15.6 


SpcX prism 


TWA 30B 


091228 


15.3 


0.69 


0.59 


1.28 


720 


± 


20 


3.1 ± 0.4 


4.4 


2.8 


19.8 


26.6 


SpcX prism 


TWA 30B 


100102 


15.1 


0.82 


1.06 


1.88 


703 


± 


8 


3.8 ± 0.2 


4.5 


2.9 


20.7 


28.5 


SpeX prism 


TWA 30B 


100103 


15.3 


0.95 


1.29 


2.24 


713 


± 


6 


2.77 ± 0.12 


4.2 


2.6 


18.5 


24.2 


SpcX prism 


TWA 30B 


100124 


15.2 


0.66 


0.50 


1.16 


630 


± 


50 


8 ± 4 


5.2 


3.6 


25.8 


39.7 


SpcX prism 


TWA 30B 


100125 


15.4 


0.68 


0.50 


1.18 


750 


± 


20 


2.6 ± 0.4 


4.1 


2.5 


18.2 


23.5 


SpcX prism 


TWA 30B 


100126 


15.2 


0.65 


0.47 


1.12 


710 


± 


30 


3.7 ± 0.7 


4.1 


2.8 


20.3 


27.7 


SpcX prism 


TWA 30B 


100127 


15.0 


0.67 


0.48 


1.15 


760 


± 


50 


3.0 ± 0.9 


4.1 


2.5 


17.7 


22.6 


SpcX prism 


TWA 30B 


100128 


15.3 


0.67 


0.46 


1.13 


780 


± 


30 


2.3 ± 0.4 


3.9 


2.3 


16.8 


20.9 


SpcX prism 


TWA 30B 


100520 


15.4 


0.78 


0.84 


1.62 


700 


± 


30 


3.5 ± 0.8 


4.5 


2.9 


20.9 


28.9 


SpcX prism 


TWA 30B 


100521 


15.6 


0.80 


0.82 


1.62 


790 


± 


30 


1.9 ± 0.3 


3.9 


2.3 


16.4 


20.1 


SpeX prism 


TWA 30B 


100522 


15.5 


0.78 


0.82 


1.60 


770 


± 


10 


2.1 ± 0.2 


4.0 


2.4 


17.3 


21.7 


SpcX prism 


TWA 30B 


100523 


15.1 


0.77 


0.81 


1.58 


760 


± 


20 


2.8 ± 0.4 


4.1 


2.5 


17.7 


22.6 


SpcX prism 


TWA 30B 


100525 


16.0 


0.88 


1.02 


1.90 


800 


± 


20 


1.5 ± 0.1 


3.8 


2.2 


16.0 


19.4 


SpeX prism 


TWA 30B 


100526 


15.8 


0.86 


1.12 


1.98 


758 


± 


12 


2.0 ± 0.2 


4.1 


2.5 


17.8 


22.8 


SpcX prism 


TWA 30B 


100527 


15.8 


0.87 


1.08 


1.95 


740 


± 


10 


2.2 ± 0.2 


4.2 


2.6 


18.7 


24.5 


SpcX prism 


TWA 30A 


091227 


9.9 


0.72 


0.33 


1.05 


















SpeX sxd 


TWA 30A 


091231 


10.1 


0.75 


0.35 


1.10 


















SpeX sxd 


TWA 30A 


100102 


9.9 


0.70 


0.32 


1.02 


















SpeX sxd 


TWA 30A 


100124 


10.0 


0.72 


0.33 


1.05 


















SpeX prism 


TWA 30A 


100125 


9.6 


0.72 


0.37 


1.09 


















SpeX sxd 


TWA 30A 


100126 


10.0 


0.66 


0.29 


0.95 


















SpeX sxd 


TWA 30A 


100520 


10.8 


1.04 


0.50 


1.54 


















SpeX sxd 


TWA 30A 


100525 


12.2 


1.53 


0.84 


2.37 


















SpeX sxd 



a All magnitudes and colors are derived using 2MASS filters for spectrophotometry on the NIR spectra taken in prism mode (as indicated) 
with the exception of the epoch marked '2MASS PSC,' which is photometry. In photometric conditions, the derived absolute spectropho- 
tometry, flux c alibrated with near by AO V stars, are accurate to ~10%, while the derived relative spectrophotometry (colors) are accurate to 
a few percent llRavner et al.ll2009T) . For a full description, see Section 3.3.1. 

b UT dates are recorded as YYMMDD. 

c Magnitudes. 

d Temperatures (in K) of the NIR excess measured by fitting BB curves to the residuals of the M4 template Gl 213 subtracted from the 
TWA 30B data. 

°The apparent area of the emitting region, in units of R0 2 . 

f Values denoted by "1" are derived from Equation [1] (the annulus of a thin disk case) while those denoted by "2" are derived from 
Equation [2] (the vertical rim/disk warp case). 

g The radial distance to the emitting black body flux arising from the disk, in units of Rq . 

h The Keplerian rotation periods for each case, in units of days. 
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Fig. 1. — Finder chart for TWA 30A (marked by a white box) and TWA 30B (marked by a red 
box). The field is centered on TWA 30A and is 5 arcmin on a side, with north up and east to the 
left. This image is a composite 3-color image made from DSS II B-, R-, and /-band images. 
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Fig. 2.— The MagE spectra of TWA 3 0B from 2009 March 4 (b lack) and 2009 May 30 (red). The 
young template RECX 5 (e Cha 5, M4; iLuhman fe; Steeghsll2004l ) is shown below these data in gray 
and provided the best spectral type match over 6000-8000 A. All spectra have been normalized at 
7500 A, and the template has been shifted down by 0.5 units for clarity. The MagE spectra have 
a resolution of R ~4100 while the template spectrum has a resolution of R ~1000; however, the 
data were boxcar smoothed down to a comparable resolution for spectral classification. 
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Fig. 3. — The radial velocity profiles of several prominent features in the two MagE spectra of 
TWA 30B, color-coded identically to that of Figure [2j For each panel, the flux has been normalized 
at —250 km s^ 1 and the star's rest velocity (+12 km s _1 ) has been subtracted out. Each feature 
is labeled in the inset. 
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Fig. 4. — Left: The SpeX prism spectra of TWA 30B at the stated epochs. All spectra have 
been normalized at 1.27 fim. The spectrophotometry and colors calculated from these data are 
listed in Table [SJ The inset shows the SXD data of TWA 30B from 2010 Jan 3 (yellow) over the 
region encompassing the NIR He I and [C I] emission features. Right: The spectrum of TWA 30B 
from 2009 May 14 (black) in comparison to the field standard Gl 51 (red); both spectra have been 
normalized at 1.27 /xm. The subtracted flux of the template (Gl 51) from the data (TWA 30B) is 
shown as a dashed gray line. We have fit a blackbody with T e jj=670 K (purple) to the subtracted 
flux. 
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Fig. 5. — Near-infrared colors of TWA 30A (red five pointed stars) and TWA 30B (green five pointed 
stars) derived from our multi-epoch SpeX data, plotted along wit h the single ep och of 2MASS data 
available for each source (s quares). The dwarf track (>M0 from [Leggettlll992l; earlier types from 
Kenvon &: Hartmannl Il995l ) is the lower solid line and the giant track (jBessell Brettl Il988l ) is 
the upper solid line, the color- color space occupie d by reddened dwarfs is outlined by two parallel 
dotted tracks, the cTTS locus (jMeyer et al.lll997l ) is shown as a single dashed line, and the vector 
at the lower right indicates a reddening vector of Ay=8 mag. Both the dwarf and giant tracks 
have been transf ormed to the CIT photometric system, which is similar to the 2MASS system 
(|Carpenterll200ll ). Each set of stars is connected by a solid straight line to delineate the sequence 



of observations. Note that the TWA 30 A track is parallel to the reddening vector while the TWA 
30B track is roughly parallel to the cTTS locus. 



-29- 



100 




600 



700 800 
BB Excess Temperature (K) 



900 



Fig. 6. — The calculated area of the emitting flux region in units of R 2 as a function of the 
blackbody (BB) temperatures calculated from the excess emission for each set of data. To calculate 
the size of the emitting region we assumed that for the scaling factor (R/d) 2 , R = 1 R© and d = 42 
pc. For further details, see Section 4.1. 
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